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Simulating Observations of Dark Matter Dominated 
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ABSTRACT 

Low Surface Brightness (LSB) galaxies are dominated by dark matter, and their rota- 
tion curves thus reflect their dark matter distribution. Recent high-resolution rotation 
curves suggest that their dark matter mass-density distributions are dominated by 
a constant-density core. This seems inconsistent with the predictions of Cold Dark 
Matter (CDM) models which produce halos with compact density cusps and steep 
mass-density profiles. However, the observationally determined mass profiles may be 
affected by non-circular motions, asymmetries and offsets between optical and dy- 
namical centres, all of which tend to lower the observed slopes. Here we determine the 
impact of each of these effects on a variety of halo models, and compare the results 
with observed mass-density profiles. Our simulations suggest that no single systematic 
effect can reconcile the data with the cuspy CDM halos. The data are best described 
by a model with a soft core with an inner power-law mass-density slope a = — 0.2±0.2. 
However, no single universal halo profile provides a completely adequate description 
of the data. 
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1 INTRODUCTION 



Low Surface Brightness (LSB) galaxies are dominated by 
dark matter. Their stellar populations are dynamically 
unimportant and the rotation curves are thus direct trac- 



& Bosma 2002; Marchesini et al. 2002) show that real ro- 
tation curves rise less steeply than predicted, and do not 
have the re quired CDM s hape. Similar conclusions have b een 
reached by palucci ( 2001 ) and 



Salucci fc Borriello ( 2001 ) for 



ersof their dark matter distributions (ie Blok & McGaugh 



1997). This is important, as in brighter galaxies it is diffi- 
cult to disentangle the dynamical contributions of the dark 
and visible matter, and unambiguously determine the dis- 
tribution of the dark matter. LSB galaxies give us a unique 
opportunity to compare observations with cosmological dark 
matter simulations. 

These numerical models, based on the Cold Dark Mat- 
ter (CDM) paradigm, make very specific predictions about 



high surface brightness disk galaxies. At small radii, the em- 
pirically determined mass distribution can be well described 
by a central kpc -sized constant-den s ity core (i.e. a ~ in 



the inner parts) (de Blok et al 



2001 



de Blok & Bosma 



2002 



the distrib ution of dark matter in galaxies ( Dubinksi & Carl- 
berg 1991). In CDM halos the mass density distribution in 
the inner parts is characterized by a steep density "cusp", 
usually described b y a power-law p(r) ~ r° with slopes va ry- 

to 



ing from a 



4 



Navarro, Frenk & White 



199£ 



1997ft 



1.5 (e.g. [Moore et al||1998[ |1999j ; |Bullock et al.|p001[ ) 



This cusp manifests itself in a steeply rising rotation curve 



ies (Moore b-994 McGaue 


h, Rubin & de Blok 


200l|; 


de Blok, 


McGaugh & Rubin 2001; 


Blais-Ouellette et al 


2001 


de Blok 



[hereafter dBMBR and dBB02 respectively] 



dBMBR and dBB02 base their conclusions on obser- 
vations of rotation curves of over 60 LSB galaxies. They 
determine the mass density profiles that give rise to the ob- 
served rotation curves and fit the inner parts with a single 
power-law. The resulting distribution of inner mass-density 
slopes shows a strong peak at a = —0.2 ± 0.2 (see Fig. 2 in 
dBMBR), but with a large spread, varying from steep and 
cuspy q = —2 profiles to flat and even moderately positive 
slopes. This was interpreted as a result of resolution effects: 
dBMBR find that the best resolved curves show the flat- 
test profiles (see Fig. 3 in dBMBR and Fig. 14 in dBB02), 
contrary to what is expected for CDM. At lower resolutions 
the CDM and the ISO models both predict the same slopes, 
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leading to an a mbiguity as to which model fits lower resolu- 
tion data best (van den Bosch & Swaters 2001). 



However, resolution is not the only effect that can influ- 
ence the measured slopes. Other effects, which dBMBR and 



dBB02. These rotation curves are measured from long-slit 
Ha spectra, mostly sampling the kinematics in the inner 
parts, in many cases combined with lower-resolution Hi ro- 



tation curves samplin g the outer parts (de Blok, McGaugh 



dBB02 did not thoroughly investigate, include non-circular 



van der Hulst 



1996) 



motions, offsets between dynamical and optical centres, lop- 
sidedness, asymmetries, and slit-width. These effects all have 
the potential to decrease the measured inner mass-density 
slopes. It is thus important to investigate if the observed 
shallow slopes truly reflect the underlying mass distribu- 
tion, or whether they are determined by systematic effects. 
For example, small offsets between the dynamical centre and 
the optic al centre are sometimes observed in l ate- type barred 
galaxie s (|Pisano 1 Wilcots fc Elmegrccn| 1998 ; do VaucoulcuK 
& Free fran|[l972D . If this were a systematic property of all 
LSB or dwarf galaxies then observations aligned on the opti- 
cal centre would not necessarily probe the dynamical centre. 
Even if the galaxy in question had a steep CDM density cusp 
at the dynamical centre, the derived slope would be shal- 
lower, giving the impression of a core-dominated halo. Hence 
these systematic effects should be understood, in order to 
judge whether LSB galaxies have cuspy or core-dominated 
halos. Furthermore, once the impact (if any) of these system- 
atic effects is understood, one could hope to determine how 
much they contribute to the scatter in the observed slopes, 
and whether the range in observed slopes can be attributed 
to "cosmic scatter" or whether a single halo profile for all 
galaxies suffices. 

In this paper we investigate the importance of these 
systematic effects by deriving a large ensemble of rotation 
curves of both cuspy CDM and core-dominated (simulated) 
halos, using observed curves as a constraint. We apply the 
relevant systematic effects to these simulated curves, and 
then determine the inner mass-density slopes. Our main con- 
clusion is that none of the systematic effects investigated 
can convincingly wipe out the signature of a CDM halo. In 
other words, if LSB galaxies had CDM halos, they would 
be unambiguously detectable, even in the presence of sys- 
tematic effects. The most likely interpretation is that the 
observed core-dominated mass-distribution reflects the true 
dark matter distribution in LSB galaxies. 

The organization of this paper is as follows. Section 2 
discusses some aspects of the observational comparison sam- 
ple. Section 3 describes the simulations. Section 4 contains 
the results of simulations that include various systematic ef- 
fects. In Sect. 5 we describe some tests performed with real 
data, while in Sect. 6 we attempt to use the simulations to 
improve on the halo model used. In Sect. 7 we present our 
conclusions. 



2.1 Defining a restricted sample 

These data span a large range in resolution, inclination, and 
general quality of the rotation curves. It is therefore worth 
asking whether the data themselves may in any way be bi- 
ased by systematic effects. For example, if all the core-like 
profiles found turned out to have been derived from low- 
resolution, low-quality rotation curves, one would have valid 
reasons to question the usefulness of those data and conclu- 
sions based on them. 



McGaugh, Barker & de Blok (2003) present such an 
analysis where they analyse only those rotation curves that 
resemble NFW CDM rotation curves best, and conclude that 
there is no systematic difference in the distributions of c for 
their restricted and total samples. 

Here we make a similar comparison, but concentrate on 
the properties relevant for our comparison with simulations, 
namely distributions of c and the inner slope a. Starting 
with the complete sample from dBMBR and dBB02 (see 
Fig. we apply the following cuts: firstly, we remove all 
galaxies with inclinations i < 30° and i > 85°. This leaves 
44 galaxies. The resulting distribution is insensitive to the 
choices for the upper and lower inclination limit: we have 
tried cuts with lower limits of 25° and 40°, and upper lim- 
its of 75° and 80°, and find almost identical distributions. 
Secondly, we remove all galaxies with low-quality rotation 
curves, i.e. with a small number of independent data points, 
with large error-bars, and large asymmetries. Galaxies where 
the minimum-disk assumption is clearly not valid were re- 
moved as well (e.g. UGC 6614). The resulting sample con- 
tains 39 galaxies. Finally, as the difference between cusp 
and core is most clearly visible in the innermost part of the 
galaxies, we demand that the inner part is resolved: we re- 
tain only those galaxies that have at least 2 independent 
data points in the inner 1 kpc. This leaves us with a final 
sample of 19 galaxies, listed in Table [ljand shown in Fig. |l| 

A comparison of the different stages of pruning shows 
no systematic differences between the distributions of the 
initial full sample and final restricted sample. If anything, 
the peak of the distribution of inner slopes a has shifted 
towards a more positive value. The tail towards negative 
slopes has disappeared. The peak in the c-distribution still 
occurs near c ~ 6, and the peak at c ~ is still present in 
the final sample. 



2 DATA 

The model rotation curves we derive will be compared with 
the results presented in dBMBR and dBB02. The focus will 
be on the histograms of the inner slope a (dBMBR, their 
Fig. 2) and the diagrams of inner slope a versus resolution 
r in of the rotation curve (dBMBR, their Fig. 3 and dBB02, 
their Fig. 14) . These diagrams are based on power-law fits to 



~ 60 high-resolution rotation curves from McGaugh, Rubin 


k de Blok 


( 


2001 


), ie Blok, McGaugh & Rubin (2001) and 



2.2 Comparison with literature 

Recently. 



Swaters et al. 



2002 [hereafter SMBB] have made 



a similar analysis of a sample of 15 dwarf and LSB galax- 
ies. They reach the conclusion that the inner mass density 
slopes in their sample cover the full range — l<a<0, but 
claim that though their data prefer a shallow a ~ slope, 
they cannot rule out steep a — — 1 slopes for their sam- 
ple. The analysis presented in the previous section and in 
dBMBR suggests that a large fraction of the steep slopes 
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one finds are actually due to insufficient resolution. It would 
thus be interesting to see if these effects are also present 
in the SMBB data. We subjected their data to the same 
quality criteria that we used for ours. In practice this meant 
that 5 galaxies with less than 2 independent points in the 
inner kpc were rejected (F563-V2, F568-1, F568-3, F568-VI, 
F574-1), as well as one galaxy (U572I) where the minimum 
disk assumptio n is invalid. A high-r esolution optical velocity 
field of U5721 (Garrido et al. 2002) clearly shows the pres- 



ence of a kinematical minor axis not perpendicular to the 
majo r axis. This is a classical indication of an oval distor- 
tion ( Bosma 197J). For this galaxy, there are images in the 
HST archive that show the presence of a weak bar (see also 
dBB02, their Sec. 9.2.2., noting that U5721 is also known as 
N3274). 

We are thus left with 9 well-resolved rotation curves. 
Our restricted sample has three galaxies in common with 
this restricted SMBB sample. We find that for two of these 
galaxies the measured slopes are consistent at better than la 
(UGC 731 and UGC 4325). For one galaxy (UGC 11557) the 
slopes differ significantly: SMBB find a = -0.84±0.27, while 
dBMBR find a = -0.08 ±0.23. An inspection of the density 
profiles in Fig. 1 of dBMBR and Fig. 5 in SMBB shows that 
this difference in slope is entirely due to a different choice 
of break radius, where SMBB have chosen a larger radius 
than dBMBR, resulting in a much steeper slope. Choosing 
a smaller break radius will give a value entirely consistent 
with that in dBMBR. An inspection of the break radii of 
the 8 other high-resolution SMBB profiles shows that UGC 
11557 is the only case where we would claim that the choice 
in break radius is ambiguous. We have overplotted the dis- 
tribution of slopes of the restricted SMBB sample in Fig. [lj 
The number of galaxies is small, so it is difficult to say any- 
thing about the intrinsic distribution of slopes, but it is clear 
that the SMBB sample is entirely consistent with our larger 
restricted sample (and even more so when the ambiguous 
slope for UGC 11577 is removed). 

Contrary to the conclusions reached by SMBB we do 
not find that edge-on or barred galaxies skew the distribu- 
tion of slopes. A comparison between our full and restricted 
samples shows the inclusion or omission of edge-on galax- 
ies does not shift the average value of the slope significantly 
from a ~ —0.2. We also note th at edge-on late-type galaxies 



are expe cted t o be transparent ( Bosma et al 



1992 



Matthews 



WooflEOQlD; 

furthermore the Ha spectra from which the 



rotation curves are derived do not show the broad velocity 
wings expected when projection effects play a large role. 

Conclusions regarding the importance of bars are more 
ambiguous as these depend on what one defines as a bar. 
Certainly none of the galaxies investigated here contain the 
dominant straight bars found in earlier-type galaxies. Here 
we are mostly dealing with Magellanic-type oval compo- 
nents, and it is obviously a matter of judgement to decide 
whether one deals with an intrinsically oval component, or 
a projected more circular one. 

SMBB list which of the galaxies in their sample they 
consider barred. In order to make a consistent barred/non- 
barred division of the galaxies in our restricted sample, we 
use their results to gauge and calibrate our classification. 
That any such classification remains uncertain is illustrated 
by the fact that none of the galaxies in the SMBB sam- 
ple are classified in NED as unambiguously barred (Hubble 



type SB), except UGC 2259 [Hubble type SB(s)dm]. SMBB 
classify this galaxy as unbarred, even though the bar with 
perpendicular spiral arms is easily visible on the Digital Sky 
Survey image. The inner mass-density slope found for this 
galaxy is steep at a = —0.86 ± 0.18. As bars are disk dy- 



namical features (e.g., Athanassoula 2002), they imply a sig- 



nificant disk mass. Such heavy disks more than compensate 
for any underestimate of the slope which might be caused 
by the noncircular motions they induce. The minimum disk 
assumption fails, and with it the steep halo slopes inferred 
in this extreme limit. 

Despite this UGC 2259 discrepancy, we have used the 
SMBB results to gauge which of the galaxies in our restricted 
sample should be classified as barred. We list our classifica- 
tion in Table 1; we find 13 non-barred galaxies, 5 barred 
galaxies and 1 ambigious case. 

The galaxies in our restricted sample have been selected 
not to suffer from inclination effects or projection effects. 
The majority of them is unbarred, and it is not obvious at 
all that the presence of a bar influences any of the results. 
This comparison of restricted samples shows that once a 
proper pruning is made to only include minimum-disk, well- 
resolved galaxies, most of the differences are merely a mat- 
ter of semantics. The observed shallow inner mass-density 
slopes of LSB and dwarf galaxies do therefore reflect the 
intrinsic properties of their potentials. 



3 MODELS 

3.1 Model rotation curves 

In order to model the various effects that can affect rotation 
curves we use the full expression for the observed line-of- 
sight velocity at any position in a galaxy velocity field. This 
is given by (see Fig. 0): 

V(x,y) = V + V C {R) sin i cos 9 + V cxp {R) sin i sin 6 (1) 

where V(x,y) is the observed velocity at position (x,y), 
where we define the dynamical centre of the galaxy to be 
at (0, 0); Vo is the systemic velocity, which is not relevant in 
this analysis and which will be ignored; Vc(R) is the actual 
rotation curve velocity at radius R, where 7? is measured 
in the plane of the galaxy. V cxp (-R) is an expansion velocity 
describing non-circular motions (which can be used to de- 
scribe e.g. streaming motions). Position (x,y) is measured 
in the plane of the sky. We additionally define the radius r 
as the projection of R on t he plane of the sky. As usual, i 
is the inclination. Following Begeman (1987) we define a co- 



ordinate system with the +X direction towards decreasing 
right ascension, the +Y direction towards the north, and the 
dynamical centre as origin. Position angles in the plane of 
the sky are measured from N towards E. The major axis of 
the galaxy has sky position angle <j>. Position (x, y) then has 
position angle 9 with respect to the major axis as measured 
in the plane of the galaxy given by: 

— x sin (j> + y cos (j> 



R 



y sin t 



(2) 
it cos I 

We now introduce the angle (/>' which is the sky position 
angle of position (x,y). We can express x and y as 
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x — —r cos (t^ — 4>') ~ ~ r sm < 



r sm 



r cos i 



(3) 



Combining this with Eqs. (2) gives 
(sin <j>' sin <j> + cos cj}' cos cj}) 



R 



R cos i 



(sin (j)' cos 4> — cos <f> sin < 



(4) 



Without loss of generalisation, we can make the sim- 
plifying assumption that the major axis coincides with the 
X-axis, so that (j> — tv/2, as sketched in Fig. ^| Equa- 
tion (4) then reduces to cos 6 — (r/R) sin <j>' and sin 6 = 
— (r/(R cos i)) sin <f>'. 

We are now in a position to derive the line-of-sight ve- 
locity at any position (x,y), given any rotation curve Vc and 
non-circular motion component V C x P . For example, Fig. ^ 
shows a slit positioned parallel to the major axis, but offset 
by an amount y. Assuming that V cxp = everywhere, the 
observed line-of-sight velocities along the slit are derived by 
varying a; at a constant y, for a given input Vc(R) and known 
i. The radii R = \J x 2 + (y/ cos i) 2 and r = \J x 2 + y 2 can 
be calculated, as can the position angle <j>' — arctan(t//a;), 
so that 6 can be derived. This then trivially gives the line- 
of-sight velocities. One last step we need to take, in order to 
"simulate" our ignorance regarding the systematic effects in 
this particular example, is to associate the positions along 



the slit x (rather than the radii r) with the derived line- 
of-sight velocities, because, as witn real data, these curves 
will be analysed under the assumption that no systematic 
effects are present. Other situations, including non-circular 
motions, can be modeled in a similar way, as shown below. 



3.2 Halo models 

For Vc(R) we will initially investigate two halo models, ap- 
ply identical systematic effects to both, and determine the 
resulting rotation curves and mass-density profiles. The first 
model is the so-called NFW halo, favoured by the CDM sim- 
ulations, but seemingly not by the observations. The other 
model is the pseudo-isothermal halo, preferred by observa- 
tions, but without any basis in cosmology. The properties of 
these models are summarised below. 



3.2.1 NFW halos 

The internal structure of DM halos formed in a CDM uni- 



verse wasjnvestigated in detail by Navarro, Frenk & White 



(1996, 1997). They found that the halos modeled in their 
N-body simulations could be well described by 



p(R) 



Pi 



{R/R s )(l + R/RsY 



(5) 



where R s is the characteristic radius of the halo and pi is 
related to the density of the universe at the time of collapse. 
The density and radius are strongly correlated with the mass 
of the halo. At small radii the density diverges as r _1 . This 
mass distribution gives rise to a halo rotation curve 



Vc(R) = V200 



ln(l + cx) — cx/(l + cx) 



-, 1/2 



ln(l + c) - c/(l + c) 



(6) 



where x = R/R^oo. The radius -R200 is the radius where 
t he density contrast exceeds 20 0, roughly the virial radius 
(Navarro, Frenk & White 1996). The characteristic velocity 
V200 of the halo is defined in the same way as i?200- The 
curve is characterised by two parameters: a concentration c 
and a velocity V^oo- These are not independent, but related 
through the assumed cosmology. For a standard ACDM this 
relation can be well described as 



logc= 1.191 - 0.064 log V200 - 0.032 log 2 V200 



(7) 



(cf. Navarro, Frenk & White 1997). This relationship has 



an associated scatter, but estimates in the literature dif- 
fer. Here we ad opt a conservative logarithmic scatter of 
er(logc) = 0. 18 (|Bullock ct al.||200l|) . Note that in a more 



Wechsler et al 



2002) advocate a scatter of 



recent paper 

(j(logc) = 0.14. As the precise value of the scatter is not 
critical for our results, we use the larger value in order to 
giv e the maximum amount of leeway t o the CDM models. 
See |de Blok, McGaugh fc Rubiij ( |200 1| ) for more details. 

Recent N-body simulations have found even steepe r in- 
ner mass-density slopes of a = —1.5 (Moore et al. 1998). We 
will not consider these steeper slopes here; the NFW model 
can be regarded as a limiting case. If the data do not ad- 
mit a = — 1 slopes, they certainly will not allow a = —1.5 
values. 

In Fig. ^| we show the inner part of the velocity field of 
a massless disk in a NFW halo with c = 8.6 (i. e. approxi- 
mately the typical value predicted by ACDM, see McGaugh, 
Barker fc de Blok| pOO^ ) and V200 = 100 kms _ \ The iso- 
velocity contours show a characteristic "pinch" in the very 
inner parts, which is a signature of the NFW profile. 



3.2.2 ISO halos 

The spherical pseudo-isothermal (ISO) halo differs from the 
NFW halo in that it has a central region with a constant 
density. Its density profile is 

Po 



p(R) = 



(8) 



1 + {R/Rc) 2 ' 

where po is the central density of the halo, and Rc the core 
radius of the halo. The corresponding rotation curve is given 
by 

V C (R) = v / 47rGpoi? 2 7 [l - (Rc/R) arctan(i?/i? c )]. (9) 

The ISO halo is characterised by any 2 out of 3 param- 
eters: an asymptotic velocity Voo, a core radius Rc and a 
central density po- These parameters are related through 



Voo = y/47TG Po R 2 c . (10) 

Fig. ^| shows the inner part of a velocity field of a massless 
disk embedded in an ISO halo with Voo = 100 kms" 1 and 
Rc = 1 kpc. The iso-velocity contours are much rounder, 
and do not show the "pinch" of the NFW velocity field. 

3.3 Simulations 

In order to be able to directly compare the simulations with 
the data it is important that the simulated halos are "real- 
istic" , in the sense that the range of halo parameters for the 
NFW model needs to be similar to that found in the simula- 
tions (which are constructed to resemble the real Universe), 
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whereas the ISO halo parameters need to resemble those de- 
termined from observations. We therefore use the following 
boundary conditions in the models. 

(i) halo parameters. For the NFW halos we take a ran- 
dom value for V200 between between 20 and 490 kms -1 , 



and compute the corresponding c- value using Eq. 7, mod- 



this did not change any of the conclusions appreciably. We 
assume a spherical halo and minimum disk, i.e. the dynami- 
cal significance of the stellar disk is assumed to be negligible 
and the stellar mass-to-light is set to zero. For LSB galax- 
ies minimum disk is in general a good assumption (e.g. ^ 



Blok & McGaugh 



1997). It is difficult to estimate what ef- 



ified by a random scatter as described above. To simulate 
ISO halos we choose a random value for Voo between 30 and 
300 kms -1 , roughly corresponding to the observed range 
in galaxy maximum rotation velocities. For po we choose a 
random val ue from the uniform range log(po/[A^o P c ~ 3 ]) = 
— 1.7±0.4 ( Firmani et al. 2000). The core radius Rc is then 
determined using Eq. |lo], 

(ii) inclination. A random axis-ratio is chosen in the 
range from 0.17 to 0.76, corresponding to inclinations be- 
tween 40° and 80°. The inclination has no direct effect on 
the slopes, but serves as a scale factor in the error-bars de- 
termined below, and hence the uncertainty in the slopes. 

(iii) physical resolution. For each simulated galaxy the 
resolution in kpc at which it will be "observed" is determined 
as follows. Figure ^| shows the logarithmic distribution of ob- 
served resolutions (radii r-j n of the innermost points of the 
rotation curves) from dBMBR and dBB02. This distribu- 
tion is well approximated by a Gaussian with an average 
/i(logrin) = —0.45 and a dispersion cr(logri n ) = 0.41. The 
2a spread in resolution thus varies from 0.05 kpc to 2.3 kpc, 
reflecting both the spread in distance as well as the mix of 
optical and Hi curves used. Though the precise parameters 
of this distribution are not critical, choosing a resolution dis- 
tribution that is approximately equal to the empirical one 
means we can make a direct comparison between the simu- 
lations and observations. 

(iv) error-bars. To simulate the uncertainties in the 
data, we de-correct the rotation curve for inclination and 
add error-bars. Figure ^ shows the distribution of the ob- 
servational uncertainties in the dBB02 data, not corrected 
for inclination. The best-fitting Gaussian has an average 
/u(logAF) = 0.62 9 and a dispersion o-(logAV) = 0. 437. 
As in dBB02 and |de Blok, McGaugh fc Rubin| (|200l[) we 



impose a minimum error of 4 kms _± on the data. We draw 
error-bars from this distribution and assign them to the data 
points in the simulated rotation curve. The curve is then 
corrected back for inclination, thus enlarging the error-bars 
and uncertainty accordingly. As the intensity of the Ha line 
is generally inversely correlated with the magnitude of the 
error-bar, this procedure assumes a random distribution of 
Ha intensities. This is likely to be simplistic, but any im- 
provement involves assuming a specific Ha intensity distri- 
bution, meaning the models would loose some of their gen- 
erality. 

(v) slit width. Equation [l] gives the line-of-sight veloci- 
ties along an infinitely thin slit. In reality, the observations 
have used a finite slit with a width of ~ 1". We therefore 
integrate Eq. |l| over a width of 1" and determine the average 
value of V(x,y) at each position along the thick slit. 

(vi) mass density profiles. The mass density profiles are 
determined exactly as described in dBMBR, with one dif- 
ference: dBMBR choose a "break radius" to determine the 
inner slope. Here we determine the slope by performing a 
weighted fit to the innermost 3 points. The uncertainty in 
the slope is determined in the same manner as in dBMBR. 
We have tested fitting to the innermost 2, 4 or 5 points but 



feet the stellar disk has on the halo, but most likely the halo 
becomes more concentrated under the influence of the disk. 
In the absence of applicable models we will not explicitly 
simulate this, except to note that because of the dynamical 
insignificance of the disk the effect is likely to be small. The 
minimum disk assumption also gives a hard upper limit on 
the value of the slope, as taking into account any contri- 
bution by the stars to the mass-density slope reduces the 
implied dark matter slope, independent of whether any con- 
traction has taken place or not. 

(vii) small-scale deviations. The model rotation curves 
can by necessity be nothing more than pale imitations of 
those of real galaxies. One of the features of real rotation 
curves are the "bumps and wiggles" , caused by, among oth- 
ers, streaming motions along spiral arms, random motions of 
H 11 regions, and gas motions due to the effects of star forma- 
tion. The magnitude and importance of these effects are of 
course different for different galaxies. They do undoubtedly 
play a role in the observed LSB galaxies rotation curves, but 
without modeling individual galaxies in detail it is impossi- 
ble to determine and apply these effects to the models in an 
unbiase d an d general manner. We will return to this matter 
in Sect. 4.2. It should be noted though that the many ways 



in which rotation curves can be derived are all designed to 
remove the effects of these small-scale deviations: for exam- 
ple, the tilted ring procedure averages velocities over rings 
at a certain radius, whereas long-slit rotation curves use the 
independent information from radii at the approaching and 
receding sides of the galaxy to derive the underlying global 
rotation velocities. Our omission of these deviations means 
that the uncertainties we derive in our simulations will be 
smaller than found in the data, but this is a small price 
to pay not having to introduce ill-constrained random mo- 
tions that may bear no resemblance with reality. A specific 
example will be discussed in Sect. 4.2. 



4 RESULTS 

4.1 Comparison with models 

For each set of input parameters described below we com- 
pute 600 rotation curves per halo model. This is significantly 
more than the 61 observational data-points available, but 
avoids small number statistics affecting the simulated his- 
togram. To compare the data and the simulations we scale 
the maximum of the data histogram to that of the simu- 
lation histogram. To indicate the uncertainties in the data 
histogram we over-plot an error-bar on each observed his- 
togram bin which reflects the uncertainty due to counting 
statistics in the observational sample (these error-bars are 
also scaled) (see Fig. |^ for examples). In the diagram of 
inner slope a versus resolution ri n we show observational 
data points in the background, with the simulation results 
over-plotted. To avoid cluttering the plot we do not always 
show all 600 galaxies, but select 61 galaxies (i.e. equal to the 
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number of observed data points) at random (see Fig. ^| for 
examples) . 

In order to isolate the investigated systematic effects 
from resolution effects, we distinguish in the histograms be- 
tween well-resolved galaxies (2 or more independent points 
per kpc; cross-hatched histogram) and less well-resolved 
galaxies (less than 2 independent points per kpc; single- 
hatched histograms). The cross-hatched histograms can be 
directly compared with the well-resolved sample in Fig. 

4.2 The Null-case 

We first investigate the case of observations without any sys- 
tematic effects. That is, we determine the major-axis curves 
defined by Eq. [l] and the ISO and NFW halo models, and 
follow the procedure described in Sect. i.3. We assume there 



are no non-circular motions and no mismatches of any kind. 
The results are shown in Fig. ^| As expected, the derived 
slopes are close to the theoretical values. The NFW model 
shows a pronounced peak at a — — 1. The finite width of the 
slit hardly affects the measured slopes. Similar conclusions 
apply to the ISO model. The high-resolution ISO models 
show a clear a = slope, while the least-resolved models 
show steeper slopes, comparable to those found in the NFW 
model at similar resolutions. 

The histograms show only partial agreement with the 
data. The peak occurs at the wrong value for the NFW 
model, and the histogram has the wrong shape. For the ISO 
model the peak occurs at roughly the correct value, but the 
shape of the histogram differs significantly from the data. 
Clearly some sort of effect, whether cosmic or systematic, 
is needed to make the simulations more consistent with the 
data. It is clear though that this effect needs to be more 
pronounced for NFW than for ISO. 



4-2.1 Small-scale deviations 
The simulated error-bars in Fig. ^ are smaller than those 



of the data. As discussed in item (vii) in Sect. 3.3, this is 
most likely due to the fact that we do not simulate small- 
scale features in the rotation curves. This does however not 
affect the conclusions we can draw from the simulations. Af- 
ter all, when averaged over ~ 60 rotation curves, we can 
regard the "bumps and wiggles" as essentially random: they 
introduce scatter but are expected to average out and will 
not affect the shape in a systematic way. Another reason not 
to include the small-scale deviations is that every galaxy is 
unique and has features in the rotation curve that depend 
on the positions of the Hn regions, spiral arms, star forma- 
tion etc. Given the low star formation rates, small number 
of Hn regions, and lack of pronounced spiral arms, it is un- 
clear how important any of these effects are in LSB galaxies 
compared to "normal" galaxies. It is therefore difficult to 
derive a universally applicable description short of modeling 
each galaxy individually, which is clearly unfeasible. 

We illustrate the fact that small scale deviations most 
likely do not affect the conclusions, but merely introduce 
more scatter, by investigating a naive (and probably too 
simplistic) description for the deviations which we then ap- 
ply to the simulations presented above. Fig. |^ shows a his- 
togram of the deviations of the raw observed rotation curves 



with respect to the smooth ro tation curves as presented in 
Fig. 7 of dBB02 and Fig. 1 of |de Blok, McGaugh fc Rubin 
( |200 ip . The deviation histogram is well described by a func- 
tion of the form cxp — | AV|. If the deviations in the rotation 
curves were radially uncorrelated we could add deviations 
derived from this distribution to our model curves, use this 
modified model curve to determine the underlying smooth 
rotation curve and subsequently measure the mass profile. 

However, in real galaxies the small-scale deviations are 
correlated: if e.g. streaming motions along a (resolved) spiral 
arm cause a data point to lie above the average curve, then 
the data point next to it is more likely to lie above it as well. 
Furthermore, the deviations are defined with respect to the 
underlying curve, in that the latter is derived by minimising 
the former. Ideally, we would have to introduce an extra step 
in the simulations to adjust the model rotation curves to 
the effect of the deviations. On the other hand, as we would 
afterwards be fitting to the underlying curve as defined by 
a minimisation of deviations, and not to the raw deviations, 
their impact on the results would have to be significantly less 
than suggested by Fig. |^. ft is thus not clear whether taking 
small-scale deviations into account would have a significant 
effect. 



As a crude attempt to show these effects, we have ap- 
plied to the two simulations described above, small-scale de- 
viations drawn from the distribution shown in Fig. |^. In 
order to introduce a degree of correlation between the devi- 
ations, we boxcar-smooth over four data points. With real 
data the next step would have been to fit an underlying 
smooth curve by minimizing the residuals, followed by the 
derivation of the mass profile from the smooth curve. To 
determine the mass profiles we thus prefer not to use the 
raw data, since the mass profile derivation uses the gradient 
of the rotation curve, fnstead we introduce a bit of bin- 
ning and smoothing to determine the underlying gradient. 
In practice, we thus have to soften the impact of the derived 
(unbinned) deviations on the simulated curves. Here we (ar- 
bitrarily) divide them by a factor of four before adding. The 
number of data points to smooth over as well as this soft- 
ening factor were found empirically by demanding that the 
error bars in the simulated slopes match those of the data. 
Fig. ^| shows the results. The method has introduced scatter, 
but not changed the conclusions that were derived from the 
simulations presented in Fig. |(| 

It is clear that the method described above is unsatis- 
factory. The smoothing and softening factors are arbitrary 
and ill-constrained and the amount of scatter introduced de- 
pends strongly on the values one assumes for them. Again, 
the lack of a universally applicable description of what are 
essentially stochastic and random processes makes this an 
unfeasible exercise. Furthermore, as the derivation of rota- 
tion curves uses extra information such as continuity and 
symmetry between approaching and receding sides it is not 
clear whether the net impact would be significant. There- 
fore, rather than trying to fine-tune a clearly unsatisfactory 
description we prefer to proceed without it, keeping in mind 
that there is this missing ingredient, but that it is not likely 
to affect the conclusions in a significant way. 
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4.3 Lopsidedness and asymmetries 

One of the assumptions made in rotation curve analysis is 
that the velocity field is symmetrical around the dynami- 
cal centre. This is not necessarily the case for all galaxies; 
some seem to have a "kinematical lopsidedness" (see e.g. 

1999 for an extensive description), where the 



Swaters et al 



rotation curves of the approaching and receding sides each 
ha ve diff erent slopes and amplitudes (cf. Fig. 2 in Swaters et 



al. 



199£). The amplitudes of the curves can differ by up to 



20 per cent. Interestingly, kinematical lopsidedness does not 
seem to automatically imply a spatial offset between optical 
and dynamical centres. Kinematical lopsidedness can thus 
affect the interpretation of mass-density slopes. 

Here we investigate the effect of kinematical lopsided- 
ness on the rotation curves of NFW and ISO halos. We 
model the lopsidedness by computing a rotation curve, but 
choose an incorrect point of symmetry (or incorrect systemic 
velocity) . After flipping the rotation curve around this incor- 
rect point of "symmetry" we thus end up with approaching 
and receding side curves that have different amplitudes and 
shapes, and therefore different inner slopes. We then pro- 
ceed to take the average of the inner curves, and determine 
the mass-density profile. 

As the derived slopes are steepest when the point of 
symmetry is chosen correctly, this procedure will flatten the 
measured slopes. Figure ^| shows an example of a "kinemat- 
ically lopsided" NFW curve. This procedure does not take 
into account the obser ved behaviour of th e curves at large 
radii (see e.g. N4395 in swaters et al. 1999| ; here one side of 
the galaxy exhibits a rising rotation curve, while at identical 
radii on the other side it has already flattened), but as our 
interest is in the inner slope, this will not affect any of the 
conclusions. 

We assume offsets of 0.1, 0.5 and 1 kpc. Note that 
though the magnitude of the effect is expressed in kpc, we do 
not imply that this involves a shift in dynamical centre. Our 
notation is merely a convenient way to express the incorrect 
choice of systemic velocity that we use to create the kinemat- 
ical lopsidedness. The offsets used here result in average in 
velocity offsets of ~ 15, ~ 30 and ~ 50 kms -1 respectively 
(see bottom row in Fig. [lo]). We also assume that there are 
no additional offsets between optical and dynamical centres. 

Figure [uj shows the results for the NFW halos. We see 
that lopsidedness is not important as long as the effect oc- 
curs on scales smaller than the resolution. A lopsidedness of 
between ~ 0.5 and ~ 1 kpc (which translates into an av- 
erage velocity shift of ~ 40 kms -1 ) is needed in order to 
bring most of the NFW points into the a — part of the 
diagram. The peak at a = —1 in the 0.1 kpc case is caused 
mainly by well-resolved galaxies. Small offsets are thus in- 
consistent with NFW models. This peak is still present in the 
0.5 kpc case, but is now mainly caused by lower resolution 
galaxies. The well-resolved galaxies are distributed more like 
the observed histogram, but with an excess of galaxies with 
positive slopes. The spread in slopes at these higher resolu- 
tions is much larger than observed. Large offsets of ~ 1 kpc 
are inconsistent with the observed distribution as nearly all 
well-resolved galaxies exhibit too positive slopes. 

The bottom-row in Fig. |l^ show the distribution of the 
velocity differences between approaching and receding side 
(not corrected for inclination) for the various offsets. Also 



shown is the distribution of the same differences but as a 
fraction of the halo velocity V2oo- 

In order to explain the data with NFW models, one 
thus needs to assume a net lopsidedness for all galaxies of 
~ 30 to ~ 60 kms" 1 . Estimates of the fraction of galaxies 
that are kinematically lo psided reach up to 30 to 50 per 
cent (Swaters et al. 1999), but with much smaller velocity 



differences. It is thus not possible to explain the data with 
CDM models in combination with lopsidedness alone. We 
do not show the ISO halo results: ISO halos already show a 
slope that is slightly more positive than the data (Fig. ^), 
adding in extra flattening due to lopsidedness results in even 
more positive slopes. ISO halos are thus not consistent with 
a large (i.e. shift in symmetry point > 1 kpc) kinematical 
lopsidedness. 



4.4 Non-circular motions 

A frequently made assumption is that the visible matter in a 
galaxy is on circular orbits, and therefore that the measured 
line-of-sight velocities represent the true rotation velocity. If, 
however, part of the measured velocity is due to non-circular 
motions, it is no longer representative indicator of the mass 
content of a galaxy. Systematic non-circular motions are 
most commonly observed as streaming motions along bars 
or spiral arms (we do not consider local non-circular mo- 
tions due to e.g. star formation or supernovae, as explained 
in Sect. 3.S and |l.2| ). Bright galaxies can show large stream- 
30 kms -1 in M81 (|Adler fe Westpfah] 



ing m otions: e.g. 



1996) or ~ 20 kms -1 in NGC 1365 (Jorsater fc van Moorsel 



1995). These galaxies are however mostly grand-design spi- 



rals with pronounced spiral structures. Streaming motions in 
late-type spirals such as M33, NGC 300, NGC 925 or NGC 
1744 are found to be much smaller, and barely distinguish- 
able from random motions in the disk (i.e. 5.10 kms -1 ) 
^eu^&^ar^derJHidst 1987 ; |Puche, Carignan fc Bosma 199C ; 



Pisano, Wilcots & Elmegreeri 1998| ). LSB galaxies are sim- 
ilar to these late-type galaxies and would not be expected 
to have large circular motions. Hi velocity fields also show 
no evidence for large systematic non-circular motions > 10 
kms -1 (|de Blok, McGaugh fc van der Hulst 



de Blok, Walter| ' |2002j ) 



1996 



Wcldrakc 



We model non-circular streaming motions by introduc- 
ing a radial velocity component (see Eq. ^). The effects of 
streaming motions on a rotation curve are strongest when 
they are oriented parallel to the minor axis. We approximate 
this by adding a constant velocity component to the rota- 
tion curve over the inner one-third of its radius (the exact 
range is not important as we use only the inner points to 
determine the mass-density profile) . We thus use Eq. [j] but 
change the phase of the radial component by n/2. Again we 
only show the results for the NFW halos, as the ISO halos 
are only consistent with no or very small non-circular mo- 
tions, for the same reason as with kinematical lopsidedness: 
their profiles are already slightly flatter than the data, and 
more flattening is not needed. 

We have computed models for streaming velocities of 5, 
10 and 20 kms -1 . Figure [ll| shows that systematic stream- 
ing motions of order ~ 20 kms -1 are needed to at least 
partially explain the data for the well-resolved galaxies with 
a NFW halo. There is however still an excess of data points 



8 de Blok, Bosma & McGaugh 



with slope a = — 1 formed of low-resolution galaxies, as well 
as an excess of positive slop es due to well-resolved galaxies. 
The bottom row in Fig. nil shows that for a significant frac- 
tion of the galaxies a streaming motion of this magnitude 
amounts to ~ 20 per cent or more of the halo velocity Vioo- 
In order to use streaming motions in combination with NFW 
halos, one thus needs to assume that all galaxies have sys- 
tematic non-circular motions over a large part of the disk of 
order ~ 20 kms -1 . This contradicts the observational data 
for late-type galaxies. 



4.5 Mismatched position angles 

A mis-match between the true position angle and the ob- 
served position angle does not lead to a change in slope 
for axisymmetric potentials. As can be seen from Eq. [I] it 
will merely introduce a scale-factor of order (cos 6/ cos 8') 
where 6 and 9' are the true and mis-matched position an- 
gles respectively. In the absence of non-circular motions and 
in axisymmetric potentials one can even correct for a mis- 
matched position angle. If the mass-distribution in the cen- 
tral parts is not axisymmetric, then the precise value of the 
position angle becomes important. However, in these cases 
the inner parts of the galaxy are likely to be dominated by 
the stellar component, or the halo is very triaxial. In these 
cases the minimum disk and spherical halo assumptions are 
no longer valid. One needs to make self-consistent multi- 
component galaxy models, and can no longer use the simple 
one-component models we have been using. 



4.6 Offsets between dynamical and optical centres 

It is generally assumed that in a galaxy the optical centre 
and the dynamical centre coincide. If, however, this is not 
the case, then, at least for NFW halos, measuring the kine- 
matics around the optical centre will lead to the impression 
that the slope of the inner mass density profile is shallow 
rather than steep. 

It is known that in some late-type (barred) galaxies the 
optical centre can be l ocated up to ~ 1.5 kpc away from 
the dynamical centre ( | pisano, Wilcots fc Elmegreen 1998| 



de Vaucouleurs & Freeman 1972; Weldrake, de Blok, Walter 



2002). This "morphological lopsidedness" is different from 
the kinematical lopsidedness discussed before, in that as long 
as the stellar population is unimportant, the potential can 
still be dominated by the halo and be axisymmetric around 
the dynamical centre. It should be stressed though, that for 
the majority of bright galaxies where the kinematics have 
been investigated in d etail, the optica l centre and dynamical 
centre do coincide (cf. Begeman 1987). In these galaxies the 



slit would sample the true dynamical centre. 

To model the effect of offsets between dynamical and 
(observed) optical centers, we run the simulations assuming 
a distribution of offsets and observe the halos with a slit 
offset from the dynamical center (see Fig. ^). We assume 
Gaussian distributions of offsets with dispersions a = (0.5, 
1, 2, 3, 4, 5, 6) arcseconds. For each of these seven cases we 
run the simulations twice, once producing 600 NFW halos 
and once producing 600 ISO halos. As an offset between dy- 
namical and optical center is equivalent to a slit offset, we 
draw a random slit offset from the relevant Gaussian offset 



distribution for each halo (i.e. the first case produces 600 
NFW and 600 ISO halos, each observed with offset slits, 
where the distribution of slit offsets was Gaussian with a 
dispersion a — 0.5", and so on for all other values of the 
dispersion). We have also modeled uniform distributions us- 
ing widths identical to the dispersions listed above, with 
identical conclusions as the Gaussian distributions. 



4.6.1 NFW halos 

In Fig. |l^ a random selection of 10 rotation curves derived 
assuming a Gaussian offset distribution with a — 5" is 
shown in comparison with the corresponding no-offset rota- 
tion curves. Also shown are the corresponding mass density 
profiles. For the average distance of the observed sample of 
~ 73 Mpc, this dispersion corresponds to a physical disper- 
sion of ~ 1.8 kpc. We use such a large offset dispersion to 
clearly show the differences between the curves, which would 
hardly be visible for smaller dispersions. 

The ri n — a comparison is shown in Fig. [H| The ob- 
served and theoretical distributions are distinctly different, 
as also seen in the histograms in Fig. [l4|. The NFW models 
show a pronounced peak at a = — 1 not seen in the data. 
The high-resolution and low-resolution histograms look very 
similar. The only way to explain the observed distribution 
with NFW halos and centre offsets is to assume that for 
all galaxies observed the position of the dynamical centre is 
offset from the optical centre by a ~ 3 — 4". At higher disper- 
sions the distributions becomes too biased towards a — 0. 
We can put this in context by comp aring with the average 
scale length of th e F-LSB galaxies in de Blok, van der Hulst 



k, Bothun (1995), which is 11.0". Compared to the sizes 
of the optical disks, the offsets needed are thus significant. 
For uniform offset distributions similar conclusions can be 
derived. 



4.6.2 ISO halos 

In Fig. [l^ a random selection of 10 rotation curves derived 
assuming the ISO model and a Gaussian offset distribution 
with a — 2" is shown in comparison with the corresponding 
no-offset rotation curve. Also shown are the corresponding 
mass density profiles. The r- m — a plot s and histograms for 
the ISO model are presented in Fig. [16l The ISO simulations 
are a better (though not perfect) match to the data than 
the NFW simulations. ISO halos are fairly insensitive to the 
effects of moderate offsets between optical and dynamical 
centres. We thus only present the results a < 2". 

The peak in the simulated distribution occurs at a — 0, 
offset from the observed a = —0.2 ±0.2. There is an ex- 
cess of steeper slope galaxies in the observations. The inner 
parts of real galaxies clearly have density profiles that are 
not precisely flat. A possible explanation is that the mini- 
mum disk assumption is slightly inappropriate, and that the 
stellar disk causes the steeper observed slopes. However, as 
stated before, modeling this effect is beyond the scope of 
our simple one-component models. The values of the slopes 
are insensitive to the precise value of the offset, indicating 
that if galaxies do have ISO halos, they will always be obvi- 
ous from the observations. The simplest explanation is that 
offsets between optical and dynamical center are small. 
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4-6.3 Observational offsets 

A mismatch between the centres has the same effect as an 
incorrectly centered slit. The absence of systematic differ- 
ences thus also implies an absence of these observational 
offsets. This is to be expected from technical point of view: 
the optical centres in LSB galaxies are often visible in the 
slit camera, and in most cases used to line up the slit of 
the spectrograph with the galaxy. The accuracy of a typical 
telescope pointing system and the accuracy of the offset pro- 
cedure from nearby stars is generally very good. Offset and 
pointing tests performed during the observations generally 
are repeatabl e down to the ~ 0.3" level . This was already 
addressed by |McGaugh, Rubin fc de Blok| fl200l| ) and dBB02 
who present rotation curves of a few LSB galaxies that have 
been observed multiple times at different dates on different 
telescopes by different observers an d show that these agre e 
very well with each other (see also Marchesini et al. 2002). 
Residual observational uncertainties could be introduced by 
the process of measuring the positions of the centres of the 
galaxies from e.g. CCD images, catalogs or H I observations. 
These offsets can for obvious reasons only be a small frac- 
tion of the size of the galaxy, and are therefore ~ l" or less, 
consistent with the scatter implied by the simulations. 



5 TESTS WITH REAL DATA 

To test whether some of the effects described above are ob- 
servable in practice, we re-observed one of the galaxies from 
dBB02. We chose UGC 4325, which is a typical late-type 
dwarf galaxy. It is nearby (D = 10.1 Mpc), resulting in a 
linear resolution of ~ 50 pc arcsec^ 1 . Its Hi velocity field i s 
regular, symmetric and appears undisturbed ( 3 waters 199£ ) 
with an inclination of 41°. 

UGC 4325 was observed in early Feb 2002 with the 
1.93m telescope at the Observatoire de Haute Provence us- 
ing the long-slit Carelec spectrograph. The set-up, observ- 
ing procedures and data-reduction were identical to those 
described in dBB02. We obtained Ha spectra of UGC 4325 
along its major axis (PA 231°), parallel to the major axis 
but offset by +5" (0.25 kpc), parallel to the major axis offset 
by —5", as well as two spectra centered on the galaxy but 
with position angle offsets of +30° and —30°. 

The resulting raw rotation curves (not corrected for in- 
clination) are shown in the top row of Fig. [j^. Especially for 
the offset curves the positions of the "centres" along the slit 
are difficult to determine, and we determined the systemic 
velocity and central position by looking for the point that 
gave maximum symmetry between approaching and reced- 
ing sides. However, as the curves are almost linearly rising in 
the inner parts, the precise choice of the centre of symmetry 
does not affect the slopes. The symmetrized and uncorrected 
curves arc shown in the second row of Fig. |l7]. Following |d^ 



Blok, McGaugh & Rubin (2001) we have performed local 



polynomial fits to the folded curves in order to bring out 
the underlying shape. We also resampled the curves with a 
spacing of 4". The resulting curves are shown in the third 
row of Fig. Also shown in the bottom row are the mass- 
density profiles. 

The similarity between the major-axis curve and the 
two off-axis curves is more consistent with an ISO model 



and velocity field as shown in Fig. |§| The curves with po- 
sition angle offsets only show modest changes in the slope, 
consistent with U4325 having a fairly axisymmetric poten- 
tial, and negligible streaming motions. For UGC 4325 the 
systematic effects discussed in this paper are not obviously 
present and likely small enough to be unobservable. If UGC 
4325 is typical, then the data presented here show that sys- 
tematic effects discussed in Sect. 3 are unlikely to have a 
large effect on the observed rotation curves. 



6 DISCUSSION 

If systematics effects indeed play only a small role in deter- 
mining the observed mass-density profile, some of the scatter 
in the observed value of the slope could be due to intrinsic 
differences in the shapes of the mass-density profiles. This 
can be used to improve on the results for the ISO halo. Fig- 
ure ^ shows that for this model the peak in the histogram 
occurs at values of a that are slightly too large. Here we dis- 
cuss two additional observationally motivated models that 
were designed as an im p rovem ent on the ISO model. The 
([1995D halo (Sect. 



first one is the 



Burkert 



3.1 



one the halo model proposed by 
on LSB Hi rotation curves (Sect 



6.1 Burkert halos 



Kravtsov et al 



), the second 



(1998) based 



The halo model presented in Burkert ( 1995 ) is a compromise 
between the ISO model and CDM models: it has the con- 
stant density core found in ISO halos, but shows the steep 
drop-off in density at large radii found in CDM simulations. 
The density profile of a Burkert halo can be described as 



p(R) = 



Po 



(l + R/r )(l + (R/r ) 2 ) 
and the corresponding rotation curve is given by 



(11) 



V 



\R) = 2^Gp ^4|ln[(l+-)Wl + (-) 2 ]-arctan -1,(12) 
R I L ro \ ro J ro J 



where ro is the equivalent of a core radius and po is the cen- 
tral density. This particular model gives a good de scription 
of the rotation curves of dwarfs an d LSB galaxies (|Burkert 



1995 



Blais-Ouellette et al 
T 



2001 



Marchesini et al 



2002) 



However, the steep R~ drop-off in density in the outer 
parts, as implied by the above equations, has not been ob- 
served unambiguously, and therefore fits using this model 
mostly constrain the inner part of the curve. Here we char- 
acterise the Burkert halo in terms of its maximum rotation 
velocity, making the assumption that the observed rotation 
curves do not probe the declining p ~ R~ 3 part of the ro- 
tation curve. U ma x occurs at i? max = 3.25 ro, and we can 
express ro as 



ro 



0.8606ttGpo ' 



(13) 



We describe the selection of halo parameters for the simu- 
lated Burkert halos in Sec. 3. J . 
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6.2 Kravtsov et al. (1998) halos 

The distribution of observed inner mass-density slopes peaks 
at a value of a = -0.2 (dBMBR, dBB02). Before the opti- 
cal rotation curves us ed in those analyses became available, 
Kravtsov et al. (1998) fitted power-law models to the Hi ro- 
tation curves of dwarfs and LSB g alaxies (the latter from ^ 



Blok, McGaugh & van der Hulst 1996) and found that the 



inner rotation curves could be described with a power-law 
with a slope «hi = —0.2 (the agreement between the optical 
and Hi values is certainly an indication that the effects of 
beam smearing on the mjority of the Hi curves must have 
been minor). 



set of profiles of the form ( 


Kravtsov et al. 


1995; Zhao 


1996 


Blais-Ouellette et al. 


2001) 



p(R) = 



Po 



[q + (R/ro)-* ][1 + (iJ/ro)T]<0+°O/T 



(14) 



where we have changed the order and si g ns of some of the 
exponents as defined in Kravtsov et al. (1995) and Blais 



Ouellette et al. (2001) to be consistent with our notation. 
The above equation can be used to describe NFW halos with 
(q,a,/3,j) = (0,-1,3,1), ISO halos with (1,-2,2,2) while 
the Burkert halo has parameter set (1, —1, 3, 2). The param- 
eter a determines the inner slope, 7 determines the rate of 
turnover in the mass profile, while /3 gives the slope of the 
mass distribution in the outer parts of halos. The q param- 
eter indicates the pres ence (q — 1) or absenc e (g = 0) of a 

~~ (|l99Sj) find that the 



Kravtsov et al 



constant-density core. 
LSB H i rota tion curves from de Blok, McGaugh fc van der 



Hulst ( 1996 ) can be best described with the parameter set 
(0, —0.2, 3, 1.5). They note that the 7 and j3 parameters are 
difficult to constrain and choose values suggested by CDM 
simulations. Though we will use these values, we note that 
the observed mass profiles of LSB galaxies are equally well 
described by asymptotically flat rotation curves with (3 = 2. 
As we are here only interested in the inner parts of the halos, 
this choice does not affect the results. 

The corresponding rotation curve can be described by 
a model of the form 

(R/ro) 3 



V{R) = V t 



[1 + (R/ro)' 



Ug+b)/a- 



(15) 



The parameter g gives the inner slope of the rotation curve 
and is related to a by g — 1 + a/2 (note our definition 
as a as a negative number). We thus find (a,b,g) = (1.5, 
0.34, 0.9). The maximum rotation velocity is related to Vt 
by Vmax = V t (g/a,y /a (l + g/b)- < - b+3)/a . The parameter sets 
(g,a,/?,7) = (0,-0.2,3,2) and (a,b,g) = (1.5,0.34,0.9) 
thus fully describe the mass-density profiles and rotation 
curves of a halo with an a — —0.2 inner mass-density pro- 
file. These profiles will be used to optimise the description 
of the data by our simulations. 



6.3 Scaling the models 

The NFW and ISO models each have a well-defined range 
of allowed halo parameters, either derived from cosmologi- 
cal simulations or observations. These ranges are less well 
defined for the Burkert and Kravtsov halos. In order to sim- 
ulate these halos with realistic parameters we assume that 



the ISO curves modeled here are a good description of ob- 
served rotation curves. We scale the Burkert and Kravtsov 
models to show the same diversity in shapes and amplitudes. 

The left panel in Fig. |l^ compares the rotation curves of 
each of the three models for Voa = V max = Vt = 100 kms" 1 , 
and Rc = ro = 1 kpc. Using these identical numerical pa- 
rameters the ISO curve rises more slowly than the Burkert 
and Kravtsov curves. The latter has a lower maximum ve- 



locity, as for this model V ma x = 0.62 Vt (see Sec. 3.2). Sim- 
ply using identical halo parameters in the simulations thus 
makes the Burkert and Kravtsov halos consistently more 
compact and more likely to suffer from systematic resolu- 
tion effects. The differences in Fig. |l^ are due to the differ- 
ent scalings that ro and Rc imply in the various models and 
we can scale the curves in order to get a more consistent be- 
haviour. Equations [to] and [l^ show that for a given central 
density and maximum or asymptotic velocity, the ISO and 
Burkert radii scale as Rc/ro — 1/2.156. For a given cen- 
tral density and maximum velocity, one needs to assume a 
value of ro that is ~ 2.2 times larger than the corresponding 
value of Rc in order to get a Burkert curve that matches 
approximately the shape of the ISO curve. The Burkert and 
Kravtsov curves can be scaled in radius by considering their 
maximum rotation velocities. For the Burkert halo this oc- 
curs at r = 3.25 ro; for the Kravtsov halo at r = 1.91 ro. A 
scaling of 1.69 in radius will thus line up the Kravtsov curve 
with the Burkert one (which implies a total scaling of 3.64 
with respect to the ISO halo). 

The velocity scales can be matched with those of the 
ISO halos by multiplying the velocities of the Burkert halo 
with 0.86; that of the Kravtsov halo by multiplying with 
1.38. We do not normalise to the asymptotic velocity of the 
ISO halo as might perhaps be expected, as this would have 
resulted in substantially different shapes. The right panel 
in Fig. |l^ shows the scaled rotation curves. This procedure 
thus merely serves to get a sample of simulated Burkert 
and Kravtsov curves with shapes that match those shown 
by the ISO curves (and as observed in real galaxies). This 
minimises possible systematic effects in the simulations due 
to differences in shape or resolution. The scaling does not 
affect the values of the central slopes. 

In order to simulate a realistic range of parameters for 
the Burkert and Kravtsov halos we sample the same numer- 
ical ran ge of parameters as for the ISO halos (as described 
in Sect. 3.2.2| ), and multiply their radial and velocity scales 
by the factors described above. 



6.4 Results for Burkert and Kravtsov halos 

As we established that the systematic effects for core- 
dominated models can only be small, we concentrate here 
on the simulation results in the absence of these effects, ex- 
cept for the small a = 0.5" scatter in the position of the 
dynamical centre (see Sect. 4.6.2). Results for the Burkert 



halos are expected to be similar to those of the ISO model, 
as the Burkert model was designed to be a compromise be- 
tween a core model in the inner parts and a CDM p ~ R~ 3 
model in the outer parts. As is clear from Fig. the only 
real difference is that the Burkert models show a slightly 
larger number of steep slopes than the ISO model, due to 
the steeper drop-off in the outer parts. 

Fig. bol shows the results for the Kravtsov halos. The 
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peak of the histogram distribution now matches well (but 
was constructed to do so). The r- m — a distributions show 
slopes that are slightly too shallow for large values of n n . 
This is introduced by the rotation curve scaling we used and 
the points can be shifted horizontally by using a different 
scaling. However, as we are only interested in the global 
trends, and are not trying to model individual galaxies in 
detail, we retain our simple scaling relations. The scaling 
does not affect the inner slopes. 

We show the velocity field of a Kravtsov halo in Fig. |2lJ 
The halo parameters were chosen to give a rotation curve 
similar to the one used for the ISO velocity field in Fig. |j| 
The shape of the iso-velocity contours is very similar to the 
ISO velocity field, and do not show the "pinching" that the 
NFW iso-velocity contours show. 



6.5 Kravtsov halos with scatter 

The Kravtsov model fits the peak in the a histogram best, 
but was of course constructed to do so. One feature of the 
data that none of the three core-dominated models can re- 
produce are the wings towards steeper slopes found in the 
data histogram. It is easy to produce shallower slopes with 
systematic effects, but very few effects actually will produce 
steeper slopes. Furthermore, the "monolithic" nature of the 
models (for each model all halos have identical slopes) cre- 
ates pronounced peaks in the histograms. In real galaxies we 
might expect the contributions of e.g. the disk to differ from 
galaxy to galaxy, introducing a scatter in the values for the 
inner total mass-density slope. We test here the hypothesis 
that there is a small scatter in the intrinsic value of the in- 
ner slopes which is responsible for the observed wings in the 
data histograms. 

We use the Kravtsov model and add a small Gaussian 
scatter to the slope. We choose an intrinsic dispersion in the 
mass-density slope <r a = 0.2. This choice has no physical 
motivation, but is simply suggested by the width of the peak 
in the data histogram. Other more fine-tuned choices may 
very well be possible, our aim here is simply to describe 
the data in as few assumptions as possible with a simple 
one-component model. Adding a Gaussian scatter a a to the 
a-parameter in Eq. [w| implies modifying the ^-parameter 
in Eq. [[H] with a Gaussian contribution with dispersion a g = 

<7a/2. 

The results for the Kravtsov simulations with scatter 
a a = 0.2 are shown in Fig. The peak again matches 
(naturally). The high- and low-a wings of the histogram fit 
the data better than previous models. The value for the scat- 
ter can probably be fine-tuned, and the distribution of the 
scatter is most likely not Gaussian (our assumptions in a 
small number of cases lead to hollow halos). The observed 
distributions of slopes is thus consistent with halos with in- 
ner slopes a = —0.2 ± 0.2(ler) and no systematic effects. 



7 CONCLUSIONS 

We have investigated the effect of systematic effects such as 
non-circular motions, asymmetries and lopsidedness on the 
values of the inner slopes of various halo models. We find 
that no realistic combination of a single systematic effect 



and NFW models can explain the observed distribution of 
slopes. 

Core-dominated models, on the other hand, are consis- 
tent with the data in the absence of large systematic effects. 
The best match to the data is obtained by using halos with 
an inner mass-density slope a = —0.2 and an intrinsic Gaus- 
sian scatter in the slope of a a ~ 0.2. As our models assume 
minimum disk, these slopes should be interpreted as those 
of the combined mass-profiles of stars, gas and dark matter 
in a galaxy. 

It is perhaps possible to construct a combination of 
multiple systematic effects to explain the data in the CDM 
framework. However, such a combination of effects would 
have to completely wipe out the strong a = — 1 peak in the 
a histograms. The resulting models would be completely 
dominated by non-cosmological effects, and not contain any 
link between properties of dark matter halos and the visible 
galaxies that inhabit them. Yet, they would have to explain 
and obey the Tully-Fisher relation, for example, and thus 
require a tremendous amount of fine-tuning. 

The implied observational signature of CDM halos is 
strong and if present should be easily seen. 

We thus conclude that the trends observed in the n n — a 
plot are mostly resolution effects, in combination with intrin- 
sic scatter in the inner slopes. The cusp problem is certainly 
genuine; it can not plausibly be attributed to systematic 
errors in the data. LSB galaxies have shallow mass-density 
slopes implying that (i) the current CDM model, or more 
precisely the current generation of N-body models based on 
the CDM prescription, do not correctly describe structure 
at the scale of galaxies, or (ii) non-cosmological effects de- 
stroy the cusp. The second option reduces the elegance and 
predictive power of the CDM model. 
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Figure 1. Distribution of the NFW concentration parameter c (left panel) and the inner mass-density slope a (right panel) based on 
the data from dBMBR and dBB02. The various superimposed histograms show the different stages of pruning: open histogram: full 
sample from dBB02 and dBMBR; single hatched histogram: galaxies with i < 30° and > 85° have been removed; double hatched 
histogram: galaxies with low-quality rotation curves, asymmetries etc. have been removed; grey filled histogram: galaxies with less than 
2 independent data points in the inner 1 kpc are removed. The thick open histogram i n shows the restricted SMBB samp l e. The dashed 



curve in the left panel shows the expected distribution of c for a ACDM universe (e.g. McGaugh, Barker ■ 
200 lh. 
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Figure 2. Sketch of the geometry of an inclined rotating disk projected on the sky. All angles and distances are measured in the plane 
of the sky, except R and 8 which are measured in the plane of the galaxy. See text for more details. 
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Figure 3. Velocity fields of the inner parts of massless disks embedded in a NFW halo (left panel) and an ISO halo (right panel). The 
velocity field is seen under an inclination angle of 60°, and a position angle of 90°. The boxes measure 5x5 kpc. The vertical minor axis 
contour is kms^ 1 , increasing in steps of 10 km s~ x outwards. The NFW halo parameters arc c = 8.6 and V200 = 100 kms -1 , the ISO 
parameters are Rc = 1 kpc and Voo = 100 kms -1 . 




Figure 4. Distribution of the resolutions log r; n of the rotation curves in presented in dBMBR and dBB02. Over-plotted is a Gaussian 
with average /i(logr; n ) = —0.45 and a-(logr; n ) = 0.41 which was used as probability distribution for the simulations. 
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Figure 5. Distribution of the error-bars in the unsmoothed rotation curves presented in dBB02. Over-plotted is a best fitting Gaussian 
with fi = 0.629 and a = 0.437. This was used as the probability distribution for assigning error s. The vertical line indicates the value of 
4 kms -1 . Errors smaller than 4 kms -1 , were set to be 4 kms -1 , in a similar way as the data. 
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Figure 6. Comparison of the simulated NFW (left) and ISO (right) halos. The top row shows compares the simulated r- ln — a data 
points (black) with the observed distribution from dBMBR and dBB02 (grey data points). No systematic effects are assumed. The grey 
dotted lines converging on a = represent the theoretical variations in slope for ISO halos with Rc = 0.5, 1, 2 kpc. The lines converging 
on a = — 1 show the variation in slope for C/V200 = 9.8/50, 8.6/100, 6.1/500 (the slope only depends on the ratio c/V^oo)- The bottom 
row shows the histograms of a values for the two models. The double hatched histogram indicates well-resolved galaxies with r- ln < 0.5 
kpc. 
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Figure 7. Hist ogram of the deviations of the raw data points from the underlying smooth rotation curves presented in dc Blok, McGaugh 
& Rubii| (2001). The distribution is symmetrical, and can be well described by a function of the form N = 117.9 e~ ' v l/ u ' ut) 
by the superimposed curve. 
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Figure 8. Comparison of the simulated NFW (left column) and ISO (right column) halos with rotation curves modified by small-scale 
velocity deviations. The top row compares r; n — a data points with the observed distribution. The bottom row shows the histograms 
of a values for the two models. Apart from extra scatter, the addition of small-scale deviations has not changed the general trend. See 
Fig. | 
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Figure 9. Example of the construction of a lopsided rotation curve. The left panel shows a symmetric NFW halo rotation curve. The 
dotted lines indicate the true point of symmetry. The dashed lines indicate the "wrong" centre of symmetry used to construct the lopsided 
curve. The offset is 0.2 kpc in radius and 22 kms~ 1 in velocity. The right panel shows the lopsided curve (drawn line) which is the 
average of the approaching and receding sides' curves (dashed lines). 
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Figure 10. As Fig. ^ but now showing the results for NFW halos in combination with kinematical lopsidedness of 0.1 kpc (left column), 
0.5 kpc (centre column) and 1 kpc (right column). The bin at a = 0.5 also contains galaxies with a > 0.5. The bottom rows shows the 
distribution of the resulting velocity differences between approaching and receding sides, both in absolute terms, and as a fraction of 
halo rotation velocity V200 • The peak at V = is due to galaxies at large distance or low inclinations where the kinematical lopsidedness 
shift is much less than one resolution element. 
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Figure 11. As Fig. |6| but now showing the results for NFW halos in combination with streaming motions parallel to the minor axis. 
The left panel shows results for 5 kms -1 , centre for 10 kms - 1 and right panels for 20 kms -1 . The bin at a = 0.5 also contains galaxies 
with a > 0.5. The bottom row shows the distribution of the magnitude of the streaming motion as a fraction of the halo velocity V^OO- 



22 de Blok, Bosma & McGaugh 



200 



100 
50 



I 

(f) 









R„=0.34 _ 


f 

■■ . . 1 . 


y=4.77" : 



250 
200 
150 
100 
50 




'I 1 1 1 1 1 1 1 1 1 1 1 1 1 



«5 



y=4.76" 
■ l ... I ... I ... I ■ . 



250 
200 
150 
100 
50 




a= 5" 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



.i 



Ri„=0.80 
y=0.44" 
■ i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i ■ ■ 



2 4 6 8 10 02468 10 02468 10 02468 10 02468 10 



I 1 1 1 1 1 1 1 1 1 1 1 1 1 




R^O.14 

y=0.61" 
. I ... I ... I ... I . . 



200 


.i i i I i 


1 1 | 1 1 1 | 1 1 1 | 1 1 




400 


J 1 1 I 1 1 1 




150 








300 
















t * * j 


100 








200 


r § 




50 




R,„=0.24 




100 


i. ? 


R„=157 J 






y=1.10" 








y=i.56" : 


n 








n 


. . . i . . . 


. . . i . . . i . . 



400 
300 
200 
100 

300 



i i I i i i I i i i I i i i I i 



y=1.90" 
■ l ■ ■ ■ l ■ ■ ■ l ■ ■ ■ l ■ 



200 
150 
100 
50 




■ I i 1 1 1 1 i 1 1 1 1 1 1 1 1 



4 



*** 



f y=3.22" 
. i ... i ... i ... i ■ ■ 



200 



100 



I 1 | I I 1 | I I I | I I I | I I 



t 




R.=0.32 



y=4.90" 
I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ 



300 - 
200 t 
100 



■ 1 1 I i 1 1 1 1 i 1 1 1 1 1 1 1 1 i 



R,„=0.57 
y=9.89" 



02468 10 02468 10 02468 10 02468 10 02468 10 




log(R/kpc) 



Figure 12. A random selection of simulated NFW rotation curves (top panel) and corresponding density profiles (bottom panel). The 
open circles show the original major axis curves and profiles, while the filled circles with error-bars show the off-axis curves and profiles. A 
Gaussian scatter with a = 5" was used for these curves. The actual value is given in the sub-panels, along with the resolution (separation 
between points) in kpc. 
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Figure 14. Comparison of the simulated distribution of NFW halos slopes (hatched histogram; low-resolution galaxies: single-hatched; 
high-resolution galaxies: cross-hatched) with that actually observed (grey histogram). The dispersion of the Gaussian distribution for 
offsets is given above each sub-panel. 
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Figure 15. As Fig. [l^, but showing ISO curves and profiles. 



LSB halo profiles 




26 de Blok, Bosnia & McGaugh 



600 



550 - 



major axis 



maj.ax. + 5" maj.ax.-5" maj.ax. + 30° maj.ax.-30° 



£ 500 



450 



400 



100 




-50 50 
radius (orcsec) 



-50 50 
radius (orcsec) 



-50 50 
radius (arcsec) 



-50 50 
radius (arcsec) 



-50 50 
radius (arcsec) 




20 40 60 80 20 40 60 80 20 40 60 80 20 40 60 80 20 40 60 80 
radius (arcsec) radius (arcsec) radius (arcsec) radius (arcsec) radius (arcsec) 



100 

80 

in 60 
E 

i 40 
> 

20 



1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 



i 



i ... i . . . i ... i 



■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ ■ ■ | ■ ■ ■ | ■ ■ ■ | ■ ■ ■ | ■ ■ ■ ■ ■ | ■ ■ ■ | ■ ■ ■ | ■ ■ ■ | 



f 



+4 



J. , , 

i* 



.i« T A 







it' 







20 40 60 80 20 40 60 80 20 40 60 80 20 40 60 80 20 40 60 80 
radius (arcsec) radius (arcsec) radius (arcsec) radius (arcsec) radius (arcsec) 



I 1 1 1 1 I 1 1 1 1 I 

a = 0.27 ± 0.19 



3 -2 



-3 - 



a = -0.32 ± 0.06 -- a=0.16 ± 0.14 



J I I I I I I I l_ 



1A 



<x = -0.30 ± 0.12 



_l I I I I I L 



a = -0.07 ± 0.13 - 



_i i i i i i_ 



-10 1 -10 1 

log(R/kpc) log(R/kpc) 



-1 1 
log(R/kpc) 



-10 1 -10 1 

log(R/kpc) log(R/kpc) 



Figure 17. Observed rotation curves of UGC 4325, after applying spatial and position angle offsets. From left to right: major axis 
rotation curve; parallel to major axis offset by +5"; parallel to major axis offset by —5"; position angle offset of +30°; position angle 
offset of —30°. Top row shows the raw rotation curves. Second row shows folded rotation curves. Third row shows re-sampled rotation 
curves. Bottom row shows mass-density profiles derived from the re-sampled rotation curves. 
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Figure 18. Comparison of the rotation curves of the ISO, Burkert and Kravtsov models, before and after scaling. Left panel: rotation 
curves of the three models, assuming Voo = V ma x = 100 kms -1 , and Rq = j"o = 1 kpc. Full drawn curve: ISO halo; dotted curve: 
Kravtsov halo; dashed curve: Burkert halo. Right panel: the three curves after the scaling described in the text. 
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Figure 20. As Figs. Fuji and F3, but now showing the results for Kravtsov halos. The bin at a = 0.5 also contains galaxies with a > 0.5. 
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Figure 21. Velocity fields of the inner parts of masslcss disks embedded in a Kravtsov halo. The velocity field is seen under an inclination 
angle of 60°, and a position angle of 90°. The boxes measure 5x5 kpc. The vertical minor axis contour is kms -1 , increasing in steps 
of 10 kms -1 outwards. The halo parameters are ro = 3.64 kpc and V m ax = 138 kms -1 . 
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Figure 22. As Figs. [l| and |l4j but now showing the results for Kravtsov halos with a a a = 0.2 Gaussian scatter added to the value for 
the inner slope. The bin at a = 0.5 also contains galaxies with a > 0.5. 
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Table 1. Restricted sample 



Name 


Class. 


Name 


Class. 


F583-1 


U 


UGC 11648 


U 


F583-4 


U 


UGC 11748 


U 


cso-lv 1200211 


U 


DDO 64 


u 


eso-lv 2060140 


U 


DDO 185 


B 


cso-lv 3020120 


B 


DDO 189 


U 


cso-lv 4880490 


B 


NGC 1560 


U 


UGC 731 


B 


NGC 4395 


U 


UGC 3371 


U 


NGC 4455 


B 


UGC 4325 


U 


UGC 11583 


U 


UGC 11557 


D 







'B' denotes barred; 'U' unbarred; 'D' 
dubious 



